The 3D version of the finite element program FESTER

In this report, a detailed description of the 3-D version finite element pro-gram FESTER is given. This includes: 1. A brief introduction to the package FESTER; 2. Preparing an input data file for the 3D version of FESTER; 3. Principal stress and stress invariant analyses; 4. 2D joint element (surface contact) characterisation and its mathematical formulation; 5. Formulations of the 3D stress-strain analyses for both isotropic and anisotropic materials, plane of weakness and cracking criteria; 6. 3D brick elements, infinity elements and their corresponding shape and mapping functions; 7. Large-displacement formulations; 8. Modifications to the subroutines INVAR, JNTB, TMAT, MOD2 etc; 9. Numerical examples; and 10. Conclusions

A unified theory of calcium alternans in ventricular myocytes.

Intracellular calcium (Ca2+) alternans is a dynamical phenomenon in ventricular myocytes, which is linked to the genesis of lethal arrhythmias. Iterated map models of intracellular Ca2+ cycling dynamics in ventricular myocytes under periodic pacing have been developed to study the mechanisms of Ca2+ alternans. Two mechanisms of Ca2+ alternans have been demonstrated in these models: one relies mainly on fractional sarcoplasmic reticulum Ca2+ release and uptake, and the other on refractoriness and other properties of Ca2+ sparks. Each of the two mechanisms can partially explain the experimental observations, but both have their inconsistencies with the experimental results. Here we developed an iterated map model that is composed of two coupled iterated maps, which unifies the two mechanisms into a single cohesive mathematical framework. The unified theory can consistently explain the seemingly contradictory experimental observations and shows that the two mechanisms work synergistically to promote Ca2+ alternans. Predictions of the theory were examined in a physiologically-detailed spatial Ca2+ cycling model of ventricular myocytes

R-from-T as a common mechanism of arrhythmia initiation in long QT syndromes

Background: Long QT syndromes (LQTS) arise from many genetic and nongenetic causes with certain characteristic ECG features preceding polymorphic ventricular tachyarrhythmias (PVTs). However, how the many molecular causes result in these characteristic ECG patterns and how these patterns are mechanistically linked to the spontaneous initiation of PVT remain poorly understood. Methods: Anatomic human ventricle and simplified tissue models were used to investigate the mechanisms of spontaneous initiation of PVT in LQTS. Results: Spontaneous initiation of PVT was elicited by gradually ramping up I-Ca,I-L to simulate the initial phase of a sympathetic surge or by changing the heart rate, reproducing the different genotype-dependent clinical ECG features. In LQTS type 2 (LQT2) and LQTS type 3 (LQT3), T-wave alternans was observed followed by premature ventricular complexes (PVCs). Compensatory pauses occurred resulting in short-long-short sequences. As I-Ca,I-L increased further, PVT episodes occurred, always preceded by a short-long-short sequence. However, in LQTS type 1 (LQT1), once a PVC occurred, it always immediately led to an episode of PVT. Arrhythmias in LQT2 and LQT3 were bradycardia dependent, whereas those in LQT1 were not. In all 3 genotypes, PVCs always originated spontaneously from the steep repolarization gradient region and manifested on ECG as R-on-T. We call this mechanism R-from-T, to distinguish it from the classic explanation of R-on-T arrhythmogenesis in which an exogenous PVC coincidentally encounters a repolarizing region. In R-from-T, the PVC and the T wave are causally related, where steep repolarization gradients combined with enhanced I-Ca,I-L lead to PVCs emerging from the T wave. Since enhanced I-Ca,I-L was required for R-from-T to occur, suppressing window I-Ca,I-L effectively prevented arrhythmias in all 3 genotypes. Conclusions: Despite the complex molecular causes, these results suggest that R-from-T is likely a common mechanism for PVT initiation in LQTS. Targeting I-Ca,I-L properties, such as suppressing window I-Ca,I-L or preventing excessive I-Ca,I-L increase, could be an effective unified therapy for arrhythmia prevention in LQTS

Impact of energetic particle orbits on long range frequency chirping of BGK modes

Long range frequency chirping of Bernstein-Greene-Kruskal modes, whose existence is determined by the fast particles, is investigated in cases where these particles do not move freely and their motion is bounded to restricted orbits. An equilibrium oscillating potential, which creates different orbit topologies of energetic particles, is included into the bump-on-tail instability problem of a plasma wave. With respect to fast particles dynamics, the extended model captures the range of particles motion (trapped/passing) with energy and thus represents a more realistic 1D picture of the long range sweeping events observed for weakly damped modes, e.g. global Alfven eigenmodes, in tokamaks. The Poisson equation is solved numerically along with bounce averaging the Vlasov equation in the adiabatic regime. We demonstrate that the shape and the saturation amplitude of the nonlinear mode structure depends not only on the amount of deviation from the initial eigenfrequency but also on the initial energy of the resonant electrons in the equilibrium potential. Similarly, the results reveal that the resonant electrons following different equilibrium orbits in the electrostatic potential lead to different rates of frequency evolution. As compared to the previous model [Breizman B.N. 2010 Nucl. Fusion 50 084014], it is shown that the frequency sweeps with lower rates. The additional physics included in the model enables a more complete 1D description of the range of phenomena observed in experiments.Comment: Submitted to Nuclear Fusion 25/01/201

Another Proof of the Generalized Tutte--Berge Formula for ff-Bounded Subgraphs

Given a nonnegative integer weight $f(v)$ for each vertex $v$ in a multigraph $G$, an {\it $f$-bounded subgraph} of $G$ is a multigraph $H$ contained in $G$ such that $d_H(v)\le f(v)$ for all $v\in V(G)$. Using Tutte's $f$-Factor Theorem, we give a new proof of the min-max relation for the maximum size of an $f$-bounded subgraph of $G$. When $f(v)=1$ for all $v$, the formula reduces to the classical Tutte--Berge Formula for the maximum size of a matching.Comment: 7 page

Algebraic Rreconstruction for Parallel Imaging with Radial Trajectory

A novel algebraic reconstruction method for parallel imaging with radial trajectory is proposed. Taking advantage of projection-slice theorem, the reconstruction can be performed by first carrying out a 1-D FFT and then iteratively solving a system of linear equations. Since the data in non-Cartesian coordinate are directly used, the bulk work of gridding is not required. The feasibility of this method was demonstrated by simulation.published_or_final_versio

A shielding-based decoupling technique for coil array design

A new decoupling technique which takes advantage of shielding for coil array design is proposed. A ground plane is used to reduce mutual inductance by reducing the effective coupling loop area of circuits; in addition, a "screen track" is placed between adjacent coils to minimize magnetic interferences. This technique requires no modification of circuit and doesnít cause loss of flexibility of coil geometry. Simulation and experiments show that a ground plane beneath surface coils can reduce mutual coupling by 10-20dB, and a "screen track" can further provide a 2-8dB more isolation.published_or_final_versio

Market concentration, risk-taking, and bank performance: evidence from emerging economies

This paper investigates the relationship between market concentration, risk-taking, and bank performance using a unique dataset of the BRIC banks over the period 2003-2010. We find a negative association between market concentration and performance, in support of the “quiet life” hypothesis. We also find that banks taking a lower level of risks perform better, in favour of prudential practice. Moreover, the BRICs’ banking sectors were all negatively affected by the 2007-2008 global financial crisis with China and Russia being the least and most affected, respectively. On average Chinese and Brazilian banks outperform Indian and Russian ones, indicating that China and Brazil have more favourable institutional infrastructure. These results are robust to alternative model specifications and estimation techniques. Our analysis may have important policy implications for bankers and regulators in the BRICs and other developing and transition countries

Phylogenetic networks: A tool to display character conflict and demographic history

Evolutionary trees have the assumption that evolution and phylogeny can be represented in a strictly bifurcating manner. Firmly speaking, from one ancestral taxon, two descendant taxa emerge. Nevertheless, hybridization, recombination and horizontal gene transfer is in conflict with this straightforward concept. In such cases, evolutionary lines do not only separate from each other, but have the possibility of melting again and are called reticulations. Consequently, networks can represent evolutionary events more realistically than phylogenetic trees. Networks can display alternative topologies and co-existence of ancestors and descendants, which are otherwise not obvious when a comparison is done on several single trees or a consensus tree. Therefore, networks have the ability to visualize the conflicting information in a given data set. Moreover, the distribution, frequencies and arrangement of haplotypes in populations can reveal the phylogenetic histories of the taxa, regarding predictions from the coalescent theory. This review aims to: (1) give a brief comparison between phylogenetic trees and networks, (2) provide the overall concept of the coalescent theory, (3) clarify how phylogenetic networks can be used to display conflict data and evaluate phylogenetic histories, and (4) offer a useful starting point and guide for sequence analysis, with the aim to discover population dynamics.Key words: Phylogenetic networks, reticulation, coalescent theory, population history, character conflict

Timing Features of the Accretion--driven Millisecond X-Ray Pulsar XTE J1807--294 in 2003 March Outburst

In order to probe the activity of the inner disk flow and its effect on the neutron star surface emissions, we carried out the timing analysis of the Rossi X-Ray Timing Explorer (RXTE) observations of the millisecond X-ray pulsar XTE J1807--294, focusing on its correlated behaviors in X-ray intensities, hardness ratios, pulse profiles and power density spectra. The source was observed to have a serial of broad "puny" flares on a timescale of hours to days on the top of a decaying outburst in March 2003. In the flares, the spectra are softened and the pulse profiles become more sinusoidal. The frequency of kilohertz quasi-periodic oscillation (kHz QPO) is found to be positively related to the X-ray count rate in the flares. These features observed in the flares could be due to the accreting flow inhomogeneities. It is noticed that the fractional pulse amplitude increases with the flare intensities in a range of $\sim 2$, comparable to those observed in the thermonuclear bursts of the millisecond X-ray pulsar XTE J1814--338, whereas it remains at about 6.5% in the normal state. Such a significant variation of the pulse profile in the "puny" flares may reflect the changes of physical parameters in the inner disk accretion region. Furthermore, we noticed an overall positive correlation between the kHz QPO frequency and the fractional pulse amplitude, which could be the first evidence representing that the neutron-star surface emission properties are very sensitive to the disk flow inhomogeneities. This effect should be cautiously considered in the burst oscillation studies.Comment: Accepted by ApJ, 23 pages, 7 figures, 3 table